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We study a model for a quantum Ising spin glass in two space dimensions by Monte Carlo simula- 
tions. In the disordered phase at T = 0, we find power law distributions of the local susceptibility 
and local non-linear susceptibility, which are characterized by a smoothly varying dynamical ex- 
ponent z. Over a range of the disordered phase near the quantum transition, the local non-linear 
susceptibility diverges. The local susceptibility does not diverge in the disordered phase but does 
diverge at the critical point. Approaching the critical point from the disordered phase, the limiting 
value of z seems to equal its value precisely at criticality, even though the physics of these two cases 
seems rather different 

PACS numbers; 71.10.Nr, 75.10.Jm, 75.40.Mg 



I. INTRODUCTION 

A feature of disordered systems which has no coun- 
terpart in pure systems, is that rare regions, which are 
more strongly correlated than the average, can play a si, 
nificant role. For classical magnetic systems, Griffithsl 
showed that such regions lead to a free energy which 
is a non-analytic function of the magnetic field at tem- 
peratures below the transition temperature of the pure 
system. However, for the static properties of a classi- 
cal system, the Griffiths singularities are very weak; just 
essential singularitiesB. By contrast, Griffiths singulari- 
ties are much more spectacular for quantum phase tran- 
sitions at T = 0, especially for systems with a broken 
discrete symmetry. One model where these effects can 
be worked out in great detail, and where Griffiths sin- 
gularities dominate not only the disordered phase but 
also the critical region, is ihe one-dimensional random 
transverse-field Ising modelljQ. In that model, one finds 
very broad distributions of various quantities including 
the local susceptibility and the energy gap, and a dynam- 
ical exponent, z, which is infinite at the critical point. In 
the disordered phase, the distribution of local susceptibil- 
ities is found to be a power law, in which the power can he 
relatedQ to a continuously varying dynamical exponentEl, 
which diverges at criticality. The average susceptibility 
diverges when z > 1, i.e. over a finite region of the disor- 
dered phase in the vicinity of the critical poini, a result 
first found many years ago by McCoy and Wutl. 

Many of these surprising results, such as the power law 
distribution of susceptibilities in the disordered phase, 
are expected to hold more generalljfl. However, it is not 
clear whether the average susceptibility will diverge in 
the disordered phase for dimension, d, greater than 1, or 
whether this is a special feature of d = 1 . Here we inves- 



tigate Griffiths singularities for a two-dimensional quan- 
tum Ising spin glass system by Monte Carlo simulations. 
Additional niotivation for our study comes from exper- 
imental works on a quantum spin glass system, which 
did not find the expected strong divergence in the non- 
linear susceptibility at the quantum phase traasition. By 
contrast, subsequent numerical simulationsollj did find a 
rather strong divergence, comparable with that at the 
classical spin glass transition. Hence it seems worth in- 
vestigating whether this discrepancy might be due to 
Griffiths singularities causing the non-linear susceptibil- 
ity to diverge even in the disordered phase, thus making 
the location of the transition difficult in the experiments. 
Somewhat less detailed results on the two-dimensional 
spin glass Jiave also been reported in parallel work by 
Guo et ahliil, who, additionally, performed calculations 
in three dimensions. 



II. THE MODEL 

The two-dinmrisional quantum Ising spin glass in a 
transverse fielcO is defined via the following quantum 
mechanical Hamiltonian 



(1) 



where the {erf} are Pauli spin matrices, the Jij are 
quenched random interaction strengths and F is an exter- 
nal transverse field. A system described by this Hamil- 
tonian undergoes a quantum phase transition at zero 
temperature, T = 0, from a paramagnetic (or spin liq- 
uid) phase |tp a spin glass phase for some-critical field 
strength FcB. As is described elsewhereBllj this model 
can be mapped onto an effective classical Hamiltonian 
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in two space plus one imaginary time dimensions, with 
disorder that is perfectly correlated along the imaginary 
time axis. This classical Hamiltonian is 

n^-Y^Y. J^lS^{r)S,{T) - Y S^{t)S,{t + 1), (2) 

where the Ising spins, Silr), take values ±1, i and j 
refer to the sites on an L x L spatial lattice, while t de- 
notes a time slice, r = 1, 2, . . . , L,-. The number of time 
slices, Lt, is proportional to the inverse of the tempera- 
ture, T, of the original quantum system in Eq. (^. Pe- 
riodic boundary conditions are applied in all directions. 
The model is simulated at an effective classical "temper- 
ature" T'^^, which controls the amount of order in the 
spins. Of course, T'^', is not the real temperature T, 
which is zero at the transition, but is rather a measure of 
the strength of the quantum fluctuations. Increasing T"^' 
therefore corresponds to increasingthe transverse field F, 
in the quantum Hamiltonian, Eq. (|l|). The nearest neigh- 
bor interactions Jy are independent of r, because they 
are quenched random variables, and are chosen indepen- 
dently from a Gaussian distribution with zero mean and 
standard deviation unity, i.e. 
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where [• ■ -Jav denotes an average over the disorder. Sta- 
tistical mechanics averages for a given sample will be de- 
noted by angular brackets, i.e. (•••). The interactions 
between time slices are ferromagnetic and taken to be 
non-random with strength unity. 

The phase diagram of the model is sketched in Fig. |^. 
Because we are dealing with a two-dimensional ia-ttice, 
there is no finite-temperature spin glass transitioned, and 
the region with spin glass order therefore lies entirely 
along the T = axis in the region < T'^' ^r^c'i where 
T^^ denotes the critical point. In earlier workcl we found 
that 
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FIG. 1. The phase diagram of the two-dimensional quan- 
tum Ising spin glass. The horizontal axis can be thought of 
as one is using the effective classical Hamiltonian in 

Eq. (j^ or r if one is using the original quantum Hamilto- 
nian in Eq. (|l|). There is a critical point at T^' = For 
T"^ < there is a spin glass ordered phase. For T'^' > T^^ 
there is no spin glass order but there are Griffiths singulari- 
ties. In the region r^"^' < T'^' < T^l^j the Griffiths singularities 
are sufficiently strong that the average non-linear suscepti- 
bility diverges as the (real) temperature tends to zero. For 
average non-linear susceptibility stays finite in 
the zero temperature limit. 



where hi is a local field on site i. In the imaginary time 
formalism this can be evaluated from 
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Y'^SmS^{T)) . 



(6) 



Since the divergent response function at a conventianal 
spin glass transition is the non-linear susceptibility^, it 
is also interesting to study the local non-linear suscepti- 
bility, given by 



(loc) ^ ^/^.V 



which can be determined in the simulations from 



where 



(7) 
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III. THEORY 

Griffiths singularities arise from localized regions which 
are more correlated than the average. They do not have 
large spatial correlations, but give rise to singularities be- 
cause of large correlations in imaginary time. To focus 
on these time correlations, it is simplest to study quan- 
tities that are completely local in space, i.e. are just on 
a single site. We shall be particularly interested in the 
local susceptibility. 



d_ 

dh, 



(5) 



We consider distributions of and Xni°'^'' obtained 

both by measuring at different sites in a given sample, 
and by taking many samples with different realizations 
of the disorder. 



A. Disordered Phase 

In the disordered phase the distributions of x^^°'^'' and 
Xni°'^^ will be very broad with a power law variation at 
large values. Physically this comes from regions which 
are locally ordered. The probability of having such a re- 
gion is exponentially small in its volume, V, but, when it 
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occurs, there is an exponentially large relaxation timeQ, 
because, to invert the spins in this region at some imag- 
inary time one has to insert a domain wall of size V, for 
which the Boltzmann factor is exponentially small in V, 
as is sketched in Fig. |2| The combination of an expo- 
nentially large result happening with exponentially small 
probability gives a broad distribution of results in com- 
plete analogy to the effect of the Griffitlis phase on the 
dynamics in classical random magnetslljllj. In the latter 
case the volume to surface ratio determines the resulting 
probability distribution for the logarithm of relaxation 
times. In contrast to this the extra dimension present 
in the quantum problem gives rise to a volume to vol- 
ume ratio instead (cf. Fig. which leads to a power 
law distribution of correlation lengths in the imaginary 
time direction. The power depends on the microscopic 
details and so is expected to vary smoothly throughout 
the Griffiths phase. 




FIG. 2. The strongly coupled space region (cluster) of vol- 
ume V ^ L'^ tends to order (locally) the spins along the imag- 
inary time (r) direction, indicated by the plus and minus sign 
meaning a spin orientation parallel (plus) or anti-parallel (mi- 
nus) with respect to the ground state configuration of the iso- 
lated cluster. The insertion of a domain wall costs an energy 
E ~ aV, where ct is a surface tension (note that the couplings 
in the r-direction are all ferromagnetic). This event occurs 
with a probability exp(— B/T"^'), resulting in the exponen- 
tially large (imaginary) correlation time ~ exp(<TL'*/r'^'). 



AE = El — Eq, where Eq is the ground state energy 
of the quantum system and Ei is the first excited state. 
This small gap is responsible for the large susceptibility 
because the latter is essentially proportional to the in- 
verse of the gap because the matrix elements which enter 
T^(ioc) cJo not have very large variations. Since there is no 
characteristic energy gap it is most sensible to use loga- 
rithmic variables. Hence, if the power in the distribution 
of InAE is A, say, then we have 



P(lnA£;) = AEP{AE) 
~ L'^AE^ 



= [lae^'' 



(10) 



where the last line defines the dynamical exponent, z, in 
the conventional way as the power relating a time scale 
to a length scale. Comparing the last two expressions we 
see that 



A 



(11) 



Hence the tail of the distribution of AE has the form 



In [P(ln AE)] = - In AS + const. 

z 



(12) 



Since the local susceptibility is proportional to the in- 
verse of the gap, the power law tail of its distribution 
should be given by 
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= --lnx('' 



const. 



(13) 



The non-linear susceptibility involves three integrals 
over time, whereas the linear susceptibility only involves 
one. Hence we assume that the distribution of Xni°'^'' is 
similar to that of {x^^°^^)'^ ^ which leads to the following 
power law tail in the distribution: 
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const. 



(14) 



Hence there should be a factor of 3 between the powers 

in the distributions of In x^'°'^' and In Xrd'^'' ■ We shall see 
that this prediction is confirmed by the numerics. As a 
result, one can characterize all the Griffiths singularities 
in the disordered phase by a single exponent, z. 

The average uniform susceptibility will diverge, in the 
disordered phase, at the same point as the average local 
susceptibility because spatial correlations are short range 
and so cannot contribute to a divergence. From Eq. (03) 
we see that this happens when 



We can relate the power in the distribution to a dy- 
namical exponent, defined in the disordered phase, as 
foUowsQ. The excitations which give rise to a large x'"'^°'^^ 
at T = are well localized and so we assume that 
their probability is proportional to the spatial volume, 
L'^ . These excitations have a very small energy gap. 



z> d . 



(15) 



Similarly, according to Eq. (|i4|), the average non- linear 
susceptibility will diverge when 



d 



(16) 
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One can also infer the nature of the divergence of Xni°'^'' 
and x''^°'^^ as the (real) temperature T tends to zero, 
see Fig.^l For x'"^"'^'' '^s expect that the distribution 
in Eq. (pf) will be cut off at x^^""'' ^ which gives 



-)1 



(17) 



For Xni'^'' we expect that the cutoff will be at of order 
T^^, which, together with Eq. (|lj) gives 



(18) 



The global non-linear susceptibility will diverge in the 
same way, possibly with logarithmic corrections- as oc- 
curs in the one-dimensional random ferromagnettl. 

The dynamical exponent will tend to some limit as 
the critical point is approached. It is interesting to ask 
whether this limit will be the same as the value of z 
precisely at criticality. On the face of it, there does not 
seem any reason why they should be equal, since z in the 
disordered phase is determined by rare compact clusters, 
whereas z at criticality is determined by fluctuations on 
large length scales of order of the (divergent) correlation 
length. Nonetheless, for the 1-d random ferromagnet, 
they are both equal (to infinity). For the 2-d spin glass 
we shall also find that these two values are numerically 
close, and may well be equal (though finite). 



B. The Critical Point 

We expect that the distribution of will also have 
a power law at the critical point just as it does in the 
disordered phase. To deduce the exponent, note that the 
average time dependejit-correlation function at criticality 
is given by scaling asBll3 



[{SmS^{T))]. 



where 



d+ z-2 + r] 



(19) 



(20) 



is the order parameter exponent and v is the correlation 
length exponent. 

Since the average local susceptibility is just the integral 
of this over t, it follows that the distribution of Inx^'"^-* 
must have the same power, i.e. 
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Inx 
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const. 



(21) 



at criticality. In earlier work! we found z = 1.5,77 = 0.5, 
so numerically [3/{vz) is about 2/3. 

Integrating Eq. (|l9|) over r from to T"^, one sees 
that the average susceptibility (which is the same as the 



average local susceptibility for a model with a symmet- 
ric distribution of interactions, such as that used here) 
diverges as T ^ like 

Mav ^ T'^/— 1 , (22) 

at criticality. Similarly the-iayerage local non-linear sus- 
ceptibility will diverge likeollj 



(loc), 
1 Jav 



T 



2/3/ziy-3 



(23) 



at criticality. The global non-linear susceptibility will 
have a stronger divergence at criticalitya: 
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Ijav 



rp(2P-dv)/zy-Z 



(24) 



Note that Eqs. (g2|) to (|2J) refer to the situation in which 
T'^' is set to the critical value T^', and the real tempera- 
ture tends to zero, see Fig. |l|. 



IV. RESULTS IN THE DISORDERED PHASE 
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FIG. 3. The log of the distribution of the log of the local 
susceptibility for T"^' = 3.7 for different values of L and L-r- 
There is no significant dependence on L and the data is also 
independent of Lr at small x''°'^''- Increasing Lt seems to 
simply extend the range over which the data lies on a straight 
line. The solid line is a fit to the straight line region of the 
data and has slope —3.92 which gives z = 0.51 from Eq. (p^. 

We use standard Monte Carlo methods to simulate the 
model in Eq. (H). Except where noted, 2560 realizations 
of the disorder were averaged over. The simulations were 
done on parallel computers: a Parsytec GCell024 with 
1024 nodes (T805 transputers) and a Paragon XP/SIO 
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FIG. 4. Similar to Fig. ^ but for the local non-linear 
susceptibility. The straight line has slope —1.23 which gives 
z = 0.54 from Eq. (^^, in quite good agreement with the fit to 
the data in Fig. pi Since the slope is more negative than —1, 
or equivalently z < 2/3, the average non-linear susceptibility 
does not diverge at this point. 



with 140 nodes (i860XP microprocessors). Massively 
parallel machines with many medium-sized nodes (in 
terms of memory) are ideal for the problem considered 
here: as long as one physical system fits into the RAM 
of one processor one only has to set up a farm topol- 
ogy to distribute the initial seed for the random number 
generators and to collect the results for the different re- 
alizations at the end of the simulation. Apart from that, 
no communication between processors is needed, so the 
parallclization is perfectly efficient; the gain in speed is 
directly proportional to the number of nodes. 

The distributions of lnx(i°'=) and Inxi'r^ at T'^' = 3.7 
and 3.5 (both in the disordered phase) are shown in 
Figs. |-§ There is a straight line region for large val- 
ues as expected, which is independent of L, and the only 
dependence on Lr is that the tail extends further for 
larger L^- This is not surprising since there is a cutoff 
due to the finite number of time slices at ;)(^('°'^) = Lr and 
xiii"^^ = L^- At both temperatures one sees that the val- 
ues of z obtained from ;;^('°'^) and Xn°'^'' are in reasonably 
good agreement with each other. At T'^' — 3.7, we find 
that z ~ 0.51, from the data for and z ~ 0.54 from 

data for Xn°'^^ Hence, according to Eq. (|l^), the average 
Xni°'^^ does not diverge at T'^' = 3.7 because z < 2/3. At 
T'^' = 3.5, we find z ~ 0.71 from the data for x^^°''^ and 
z ~ 0.76 from data for Xta'^'' ■ Hence, the average X'ni'^'' 
does diverge at T'^' = 3.5. 



FIG. 5. Similar to Fig. Q but for T'^' = 3.5. The straight 
line has slope —2.78 which gives z = 0.71 from Eq. ([l^). 
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FIG. 6. Similar to Fig. ^ but for the local non-linear 
susceptibility. The straight line has slope —0.87 which gives 
z — 0.76 from Eq. (^), in fair agreement with the fit to 
the data in Fig. |5| Since the slope is greater than —1, or 
equivalently z > 2/3, the average non-linear susceptibility 
does diverge at this point. 



Fig. 1^ shows the values of z at various points in the 
disordered phase. In all cases there is good agreement 
between the estimates from the data for ;^('°'^) and Xni'^'' ■ 
From this data we find that the average non-linear sus- 
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FIG. 7. 



The dynamical exponent z, obtained by fitting 



the distributions of and Xni"' ^ Eqs. (|l|) and (|l|), 

is plotted for different values of T"^'. The estimates obtained 
from data for x''""^' are shown by the triangles and the es- 
timates from the data for X^]"^^ are shown by the hexagons. 
The two are in good agreement. The dcAted vertical line indi- 
cates the critical point, obtained in RefB and the solid square 
indicates the estimate of z at the critical point. The dashed 
line is 2 = 2/3; the average non-linear susceptibility diverges 
for z larger than this, i.e. T'^' > T^^.^ ~ 3.56. The solid curve 
is just a guide to the eye. 



ceptibility diverges, in the disordered phase, for 



where 



Xnl 



3.56 . 



(25) 



(26) 



and Tjf' ~ 3.275 from earlier workcl. Note that, according 
to Fig. and Eq. (|l|), the average linear susceptibility 
does not diverge anywhere in the disordered phase. It 
appears that the value of z precisely at criticality may 
equal the value as the critical point is approached from 
the disordered phase, even though it is not clear that 
they have to be equal. 



V. RESULTS AT THE CRITICAL POINT 

From finite size scaling, the average uniform suscepti- 
bility (which is the same as the average local susceptibil- 
ity for a model with a symmetric distribution of interac- 
tions, such as that used herei,A£aries with L and Lr at 
the critical point according toolij 



[X]av = L^'X 



where 



/3 z - d + 

X = z = 

V 2 



(27) 



(28) 



In earlier workcl we found z ~ 1.5, rj ~ 0.5, v ~ 1.0 
and /3 ~ 1.0 so (z - d + 2 - r7)/2 2i 1/2. The earlier 
work concentrated on an fixed value of Lt/L^, which 
we call the "aspect ratio". Here we investigate whether 
Eq. ( |2^ ) is satisfied for a range of aspect ratios. The 
data, shown in Fig. ^, does indeed collapse well with the 
expected values of the exponents. For L^/L^ ^ 1, which 
corresponds to a large system at finite temperature, the 
dependence on L should drop out, and so, from Lt cx 
T~^, one recovers Eq. (|2^). Using the numerical values 
of the exponents gives 



Ma 



-1/3 



(29) 



Thus the average susceptibility diverges at the critical 
point, though we have seen above that it does not diverge 
in the disordered phase. 
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FIG. 8. A scaling plot of the data for the average sus- 
ceptibility at T"^' = 3.30, which is the critical point, to 
within our errors, see Eq. (^. The plot assumes the form hi 
Eq. ( p7[ ) with the exponent values deduced in earlier worUj, 
i.e. z — 1.5, P/v = 1. It is seen that the plot works well. 
The solid line, which fits the data for Lt/L^ < 0.8 has slope 
0.5/1.5 = 1/3 as expected, since the average susceptibility 
should be independent of L in this limit. The power 1/3 
gives the divergence as T — » at criticality, see Eq. p9[). 

The average local non-linear susceptibility is expected 
to vary at the critical point as 
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FIG. 9. A finite size scaling plot of the average local 
non-linear susceptibility, according to Eq. (|3o|), at T'^' — 3.30, 
very close to the critical point. The solid line, which fits the 
data for small L^/L^ has slope 2.7/1.5 = 1.8 as expected, 
since the average local non-linear susceptibility should be in- 
dependent of L in this limit. The power 1.8 gives the diver- 
gence as T ^ 0, and this value agregs well with estimates 
from the earlier estimates of exponentu, see Eq. (p2|), 
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Witjh the numerical values of the exponents found 
earliem, one has y ~ 2.5. For Lt/L^ <C 1, which cor- 
responds to a large system at finite temperature, the 
dependence on L should drop out, and so one recovers 
Eq. (p3|). Using the numerical values of the exponents 
gives 
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-5/3 



(32) 



The data, shown in Fig|^, collapses well for Lt/L^ not too 
large provided y/z ^ 1.8, which gives a divergence, 
close to the prediction in Eq. (p2[). However, the data 
collapse is not as good for larger values of Lr with z = 
1.5. Data in this region is difficult to equilibrate, which 
may be the cause of the discrepancy. It should be noted, 
though, that a better data collapse is obtained for larger 
values of z. However, the data for [x]av does not scale 
well with a significantly larger value of z. 

We show results for the distribution of ;)(-('°'^) at the 
critical point, T'^' = 3.3, in Fig. HO. Unlike the data in 



the disordered phase, shown in Figs. ^ and |^, there is 
here a significant size dependence, with the slope of the 
tail becoming less negative with increasing L. Asymp- 
totically, the slope should be given by Eq. (^l|), which 
has the value —2/3 using the exponent values obtained 
earliem. The L = 6 data has a slope of —1.7 and the 
L = 12 a slope of —0.95 so it is possible that the slope 
would tend to —2/3 for L — > oo. However, it is also possi- 
ble that the slope might be less negative then this, which 
would imply a value of z larger than 3/2. 
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FIG. 10. High precision data with large Lr for the 
distribution of the local susceptibility at the critical point, 
T'^^ = 3.30. The number of samples was 25600 for L = 6 and 
10240 for L = 8 and 12. The dashed line, which is a fit to the 
L = 6 data, has a slope of —1.7, and the solid line, which is 
a fit to the L = 12 data, has a|-slope of —0.95. Using the val- 
ues of exponents found earlieru the slope, given by Eq. (pl[), 
is expected to be about —2/3 in the thermodynamic limit, 
and it is certainly plausible that data for larger sizes would 
extrapolate to this value. 



VI. GLOBAL NON-LINEAR SUSCEPTIBILITY 

Experimentallyi one measures the global non-linear 
susceptibility, for which the local magnetization (erf) in 
(0) and the local external field hi have to be replaced 
by the mean global magnetization, L~'^J2i{'^i) ^ 
uniform field, H, respectively. For the effective classical 
model this means that one has to consider 



1 



Xnl = 



(33) 



with M = J2i T ^ii^)- earlier worm we found that this 
quantity diverges at criticality (for fixed aspect ratio) like 
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X„i L^+'^j with y given by (pi]), so for arbitrary aspect 
ratio, finite size scaling gives 



Xnl 1 



(34) 



at criticality. Here, we have looked at scaling of var- 
ious moments of the global non-linear susceptibility at 
the critical point, for a range of aspect ratios. For ex- 
ample the average is shown in Fig. 11. Since there are 
two exponents which can be adjusted, y and z, the data 
is unable to determine them both with precision. How- 
ever, in the limit Lr/L^ ^ 1, where the dependence on 
L drops out and 



[Xn 



Ij av 



j^{d-\-y)lz ^ rp-{d-\-y)/z 



(35) 



the data constrains {d + y)/z to be about 3, in agreement 
with|-jthe divergence at criticality found in earlier 
workB. Fig. ^ assumes the previously determined value 
of z, i.e. z = 1.5. 
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FIG. 11. A plot of the average global non-linear suscep- 
tibility at criticality, for a range of sizes and aspect ratios. 
The exponents jused in this fit are z — 1/5 and y + d = 4.7. 
In earlier workH, we found z ~ 1.5, y -I- d ~ 4.5 so the 
present results are consistent with these estimates. In the 
limit Lt <C , the average global non-linear susceptibility 
varies as l'^^^^^^' ^ y-(<i+H)/z giving a strong divergence of 
roughly T'"^'^ . This behavior is shown by the solid line, which 
is a fit to the data for small Lt/L^ and has slope 3.1. 

In addition we have evaluated the typical global non- 
linear susceptibility defined by 



typ 
Anl 



exp[logXn 



Ijav : 



(36) 



at the critical point and show the data in Fig. U2. As with 
the data for the average in Fig |l^, the ratio (d + y)/z is 



FIG. 12. A plot of the typical global non-linear suscepti- 
bility, defined in Eq. (^6|), at criticality, for a range of sizes 
and aspect ratios. In the limit L-r ^ L" , the typical global 
non-linear susceptibility varies as L^^!"^^^^^ ~ y-(d-i-H)/z 
ing a quite strong divergence of roughly T~^''^ . This behavior 
is shown by the solid line, which is a fit to the data for small 
Lt/ L'^ and has slope 2.4. 



more tightly constrained than either z or d+y separately. 
A good fit is obtained with {d+y)/z ~ 2.4, which leads to 
a divergence of roughly T-^ "*, not quite so strong as from 
the average non-linear susceptibility. The difference may 
well reflect corrections to scaling for the range of sizes 
that we were able to study. Note that even the typical 
value has a strong divergence with T at criticality, in 
contrast to the experimentsO. 

We have also studied the probability distribution 
P(Xni) in the disordered phase. This shows a slightly 
more complicated behavior than the probability distribu- 
tion of local quantities presented above. We find that the 
power describing the tail in the same distribution is the 
same as that of the local non-linear susceptibility. This is 
reasonable since these unusually large values come from 
correlations which are very long ranged in time, whereas 
the spatial correlation length is small and so does not 
give a significant extra effect. These spatial correlations 
do, however, cause the peak in the distribution to shift 
to larger values with increasing system size, though pre- 
sumably the peak position would eventually settle down 
to a constant for sizes greater than the correlation length. 



VII. CONCLUSIONS 

One can characterize Griffiths singularities in the dis- 
ordered phase of a quantum system undergoing a T — 
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transition with discrete broken symmetry, by a single, 
continuously varying, dynamical exponent, z. Average 
response functions may or may not diverge in part of the 
disordered phase near the critical point, depending on the 
value of z, see Eqs. ( [l5| ) and (|l6|). The numerical results, 
summarized in Fig.]7[and Eqs. (^ ) and (26), indicate 
that the average linear susceptibility does not diverge in 
the disordered phase of the 2-c? quantum Ising spin glass, 
though it diverges at the critical point. The average non- 
linear susceptibility does, however, diverge in part of the 
disordered phase. 

Numerically, as one approaches the critical point from 
the disordered phase, the value of z is close to the value 
obtained precisely at the critical point, see Fig. 0. Since 
the same result is known to feold exactly in where 
they are both equal to infinitytl, one might speculate that 
they are equal in general, though we are not aware of any 
proof of this. Presumably the detailed dependence of z 
with T^' in the disordered phase, shown in Fig. is non- 
universal. However, it is interesting to ask whether the 
answer to the question "Does the non-linear susceptibility 
diverge in the disordered phase!'' is universal or not. From 
Eq. ( p^ ) this depends on whether z > (i/3 as the critical 
point is approached. If this limit for z is precisely the 
same as the value of z at criticality, then the answer to 
the question is universal. However, as we just mentioned, 
we are not aware of any argument which shows that these 
two values of z should be equal in general. 

A related-Study has also been carried out recently by 
Guo et alO'EZi Their results for the 2-d spin glass are 
consistent with ours, and they also performed some sim- 
ulations for the 3-d spin glass. In three dimensioos, the 
classical model has a finite temperature transitior£3, so 
the spin glass phase would exist for a finite range of T, 
which shrinks to zero as T'^' — > T^' , see Fig. 0. The 
results of Guo et alO indicate that, in three dimensions, 
the range of the Griffiths phase over which the non-linear 
susceptibility diverges, i.e. the region between T^' and 
T^Jjj in Fig. |l|, is very small but apparently non-zero. It 
is interesting to speculate on whether the possible diver- 
gence of the non-linear susceptibility in part of the dis- 
ordered phase nught be related to the difference between 
the experimentsH which apparently do not find a strong 
divergence fl^ni at the quantum critical point, and the 
simulationgalld which do. 

The data presented here is consistent with our ear- 
lier resultaj in finding a dynamic exponent at the critical 
point of about 1.5. Tlus.is rather different from the situa- 
tion in one-dimensionQ'tl where z = oo, and one might ask 
whether the true dynamical exponent might not be larger 
than 1.5 in two-dimensions, and possibly infinite. While 
the data for the modest range of sizes that can be stud- 
ied by Monte Carlo simulations is consistent with a small 
value of z, we cannot completely rule out the possibility 
that this estimate would increase if one could study larger 
sizes. Unfortunately, it does not seem feasible to study 
very much larger sizes with current computer power, un- 



less a more sophisticated algorithm can be found than 
the single spin-flip approach used here. Assuming that z 
is indeed finite, the critical scaling in the two-dimensional 
quantum Ising spin glass is of a fairly conventional, but 
anisotropic, type, with z playing the role of an anisotropy 
exponent. The difference from a classical magnet with 
anisotropic scaling is that Griffiths singularities give ad- 
ditional singularities in various scaling functions in the 
limit Lt ^ L^, or equivalently TL^ <C 1. 
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